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Calculations  are  made  of  the  ionization  potential  (IP)  and  the  oxidation  potential  (Eox)  values  of  108  organic 
molecules  that  are  potential  electrolyte  additives  for  the  overcharge  protection  of  lithium-ion  batteries 
(LIBs).  The  calculated  Eox  values  are  in  close  agreement  with  the  experimental  ones,  where  the  root- 
mean-square  deviation  is  0.08  V  and  the  maximum  deviation  is  0.15  V.  The  molecules  exhibiting  high  Eox 
(>4.5  V)  show  one  of  the  following  two  features:  (1)  IP>7.70eV  or  (2)  IP<7.70eV  with  a  relatively  large 
molecule  size.  Consideration  of  bulk  solvent  effects,  in  particular  the  electrostatic  attraction  between 
solute  and  solvent,  is  crucial  in  determining  Eox.  Considering  its  accuracy  and  reliability,  the  density 
functional  calculation  is  recommended  as  a  useful  tool  for  screening  electrolyte  additives  for  LIBs. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Overcharge  events  are  an  urgent  safety  issue  in  the  lithium-ion 
battery  (LIB)  industry.  Presently,  electrolyte  additives  are  widely 
employed  for  the  overcharge  protection  of  commercial-grade  LIBs. 
Among  them,  a  redox  shuttle  (RS)  bypasses  the  surplus  cur¬ 
rent  via  a  repetitive  redox  reaction  between  the  cathode  and  the 
anode  [1-5].  The  RS  compounds  have  their  own  redox  potentials, 
selected  to  prevent  self-discharge  or  shuttling  during  charging, 
which  should  be  0.3-0.4  V  higher  than  the  end-of-charge  potential 
of  the  cathode.  Aromatic  compounds  are  another  class  of  elec¬ 
trolyte  additives  used  for  overcharge  protection.  Polymerizable 
monomers  such  as  cyclohexyl  benzene  (CHB)  and  biphenyl  (BP) 
are  electrochemically  oxidized  at  the  overcharged  cathode  to  form 
a  passive  polymeric  film  at  the  cathode] electrolyte  interface  [6-12]. 
In  addition,  non-polymerizable  aromatic  compounds  such  as  tert- 
butyl  benzene  act  as  redox  mediators,  which  facilitate  oxidative 
decomposition  of  electrolyte  solvents  that  leads  to  gas  evolution 
[6]. 

Regardless  of  the  mechanism  of  overcharge  protection,  the  elec¬ 
trolyte  additives  should  remain  inert  during  the  normal  operation 
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of  LIBs  and  become  active  only  at  overcharge  events;  this  requires 
the  additives  to  have  an  oxidation  potential  that  is  0.3-0.4  V  higher 
than  the  cathode  potential.  For  instance,  LiCo02  reaches  the  fully 
charged  state  near  4.2  V  (vs.  Li/Li+)  and  therefore  requires  additives 
with  an  oxidation  potential  higher  than  4.5  V.  Thus,  determina¬ 
tion  of  the  oxidation  potential  is  a  crucial  step  in  evaluating  novel 
electrolyte  additives. 

By  means  of  density  functional  theory  calculations,  the  present 
work  has  obtained  the  oxidation  potentials  of  108  organic 
molecules,  which  are  possible  candidates  for  LIB  additives.  Calcu¬ 
lation  accuracy  is  assessed  through  comparison  with  experimental 
data.  Understanding  the  relationship  between  the  oxidation  poten¬ 
tial  and  the  molecular  characteristics  is  considered  to  be  very  useful 
for  developing  novel  electrolyte  additives  with  specified  oxidation 
potentials. 

1.1.  Calculation  details 

Density  functional  theory  (DFT)  has  become  a  popular  method 
for  calculating  oxidation  potentials  for  a  vast  array  of  organic 
molecules  used  in  LIBs  [13-16].  The  ground-state  structures  of 
the  molecules  have  been  fully  optimized  within  Q  symmetry  by 
means  of  DFT  methods.  Spin-restricted  and  unrestricted  schemes 
have  been  employed  for  even-  and  odd-numbered  electron  sys¬ 
tems,  respectively.  All  DFT  geometry  optimizations  were  obtained 
using  DMol3  numerical-based  density-functional  computer 
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software  [17,18]  implemented  in  the  Materials  Studio  Modeling 
3.2  package  from  Accelrys,  Inc.  The  Kohn-Sham  equation  was 
calculated  with  the  gradient-corrected  BPW91  functional:  the 
1988  exchange  functional  of  Becke  [19]  with  the  correlation 
functional  of  Perdew  and  Wang  [20].  A  numerical  basis  set  of 
double-zeta  plus  polarization  (DNP)  quality  was  employed.  Utiliz¬ 
ing  the  geometries  obtained  at  the  BPW91  /DNP  level,  single-point 
energy  calculations  were  performed  at  the  B3PW91/6-31  +G(d,p) 
level  of  theory.  The  functional  includes  a  three-parameter  adia¬ 
batic  connection  exchange  term  [21],  i.e.,  a  linear  combination 
of  exact  Hartree-Fock  exchange,  Slater  exchange  [22],  and  B88 
gradient-corrected  exchange  [19]. 

This  study  employs  the  conductor-variant  polarized  continuum 
model  (CPCM)  [23],  which  places  the  solute  in  a  molecular-shaped 
cavity  imbedded  in  a  continuum  dielectric  medium.  In  the  CPCM 
models,  the  variation  of  the  free  energy  when  going  from  vac¬ 
uum  to  solution  is  composed  of  the  work  required  to  build  a 
cavity  in  the  solvent  (cavitation  energy)  together  with  the  elec¬ 
trostatic  (solute-solvent  interaction  and  solute  polarization)  and 
non-electrostatic  work  (dispersion  and  repulsion  energy).  A  dielec¬ 
tric  constant  of  31.9  is  adopted  and  is  a  weighted  average  value 
between  the  dielectric  constants  of  ethylene  carbonate  (EC:  89.2) 
and  ethyl  methyl  carbonate  (EMC:  2.9),  because  EC/EMC  =  1/2  solu¬ 
tion  is  used  as  the  solvent  in  cyclic  voltammetry  experiments.  The 
solvation  energies  obtained  from  the  CPCM  calculations  are  free 
energies.  In  the  self-energy  of  the  solute,  the  entropy  and  zero- 
point  contributions  are  omitted  due  to  molecular  modes,  assuming 
that  they  would  make  only  a  relatively  minor  contribution  to  the 
desired  oxidation  potentials  [24,25].  All  of  the  single-point  DFT  and 
CPCM  calculations  are  carried  out  by  using  the  program  package 
Gaussian03  [26]. 

2.  Experimental  details 

Battery-grade  1  M  LiPF6  in  EC/EMC  (1/2,  v/v)  was  chosen  as  the 
base  electrolyte  solution.  Using  the  Karl-Fisher  titration  method, 
the  water  content  in  the  electrolytes  was  determined  to  be  less 
than  20  ppm.  The  free-acid  content  in  the  electrolytes  was  less  than 
50  ppm. 

The  oxidation  potentials  of  the  benzene  derivatives  were 
determined  by  linear  sweep  voltammetry  (LSV).  A  Pt  disc 
(area  =  0.02  cm2)  was  used  as  the  working  electrode.  The  reference 
electrode  was  Li  foil  and  the  counter  electrode  was  Pt  wire.  The 
potential  was  swept  from  the  open-circuit  value  to  6  V.  The  oxi¬ 
dation  potentials  reported  in  this  paper  were  obtained  from  the 
current-potential  curve  by  measuring  the  intersection  between 


Table  1 

Solvent  effects  [AG  (cation) -AG  (neutral)  in  eV]  on  oxidation  potentials  for 
molecules  presented  in  Figs.  1  and  2. 


two  extrapolated  lines  of  the  residual  current  line  and  the  increas¬ 
ing  oxidation  current  line. 

The  preparation  of  electrolytes  and  all  electrochemical  experi¬ 
ments  was  performed  in  an  argon-atmosphere  glove  box,  in  which 
water  and  oxygen  concentrations  were  maintained  below  5  ppm 
and  the  temperature  was  held  at  25  ±  2  °C. 

3.  Results  and  discussion 

The  observed  oxidation  potential  (E0x)  values,  together  with 
the  calculated  ionization  potential  (IP)  and  Fox  values,  for  108 
molecules  are  listed  in  Table  1 S  (see  supplementary  data).  All  poten¬ 
tials  are  reported  vs.  Li/Li+.  The  calculated  oxidation  potentials  are 


Fig.  1.  (a)  Ionization  potential  vs.  £0x(exp.);  (b)  £ox(c ale.)  vs.  £0x(exp.)  for  all  108  molecules. 
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in  good  agreement  with  the  experimental  data,  where  the  root- 
mean-square  (RMS)  deviation  is  0.08  V  and  the  maximum  deviation 
is  0.15  V.  The  correlation  coefficients  (R2)  of  IP  and  E0x(calc.)  with 
respect  to  E0x(exp.)  are  0.84  and  0.98,  respectively.  This  observa¬ 
tion  indicates  that  consideration  of  bulk  solvent  effects  is  crucial 
for  determining  oxidation  potentials  (see  Fig.  1 )  [27].  The  dielectric 
continuum  model,  CPCM,  provides  reliable  solvent-driven  stabiliza¬ 
tion  effects,  in  spite  of  the  main  limitation  of  PCM  methods,  i.e.,  the 
total  absence  of  explicit  solvent  molecules.  These  results  imply  that 


specific  short-range  interactions  are  insignificant  between  ionic 
solutes  and  the  first-shell  of  solvent  molecules  [28].  The  observed 
agreement  between  the  calculated  and  the  experimental  oxidation 
potentials  allows  estimation  of  unknown  or  inaccessible  oxidation 
potentials  for  organic  molecules  in  such  non-aqueous  solutions. 

Fig.  2  shows  25  high-voltage  molecules,  each  with  an  £0x(exp.) 
value  greater  than  4.50  V.  It  is  worth  noting  that  all  the  IP  values  of 
the  molecules  are  larger  than  7.70  eV,  except  in  one  case  (7.33  eV, 
molecule  22).  There  are  only  eight  molecular  groups  with  the 


5.07(8.24) :  5.05 


2 

5.00(8.15):  5.00 


3 

5.07(8.50) :  4.95 


4 

4.89(8.00) :  4.93 


9 

4.91(8.16):  4.78 


10 

4.86(8.26) :  4.73 


11 

4.62(7.73) :  4.65 


12 

4.67(8.05):  4.61 


17 

4.67(8.17):  4.58 


18 

4.60(7.73) :  4.58 


19 

4.63(7.93) :  4.58 


20 

4.70(8.20) :  4.55 


Fig.  2.  Twenty-five  high-voltage  molecules  with  Eox( exp.)  >4.50  V.  £0x(calc.):£0x(exp.)  values  reported  along  with  IP  values  in  parentheses. 
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25 

4.45(7.73) :  4.50 


Fig.  2.  ( Continued ). 


condition  of  IP  >  7.70  eV  and  fox  <  4  .50  V,  as  shown  in  Fig.  3.  In  the 
case  of  molecule  61,  for  example,  the  Eox  value  is  merely  4.16  V, 
despite  having  an  IP  value  of  7.74  eV.  The  solvation  energies  have 
been  examined  in  detail  to  understand  such  discrepancies  between 
IP  and  Fox  values  (see  Table  1 ). 

Some  remarks  need  to  be  made  regarding  the  definition  and 
physical  meaning  of  various  solvation  energies  reported  in  this 
study  [29,30].  Solvation  free  energies  are  calculated  as  the  sum  of 
four  different  terms,  namely:  cavitation,  electrostatic,  dispersion 
and  repulsion  energies 

AGsoi  =  AGcav  +  AGjnt  =  AG  cav  +  AGeiec  +  ^Gnon-elec 

=  AGcav  +  AGeiec  +  AGdiSp  +  AGrep  (1) 

In  this  expression,  AGcav  is  the  Gibbs  free  energy  of  cavity  forma¬ 
tion  in  the  liquid  solvent;  whereas,  the  following  three  terms  are 
the  electrostatic  AGeieo  dispersion  AGdisp,  and  repulsion  AGrep 


contributions  to  the  energy  change  AGint  due  to  turning  on  the 
solute-solvent  interactions  when  the  solute  is  in  the  cavity.  The 
energy  partition  scheme  employed  here  can  be  useful  to  analyze 
the  nature  of  the  solute-solvent  interaction.  The  solute-solvent 
electrostatic  term  is  found  to  be  the  main  interaction  between 
the  solute  and  solvent  molecules,  as  shown  in  Table  1.  The  cav¬ 
itation,  dispersion,  and  repulsion  contributions  to  the  oxidation 
potentials  play  only  a  minor  role,  because  such  effects  for  cation 
and  neutral  molecules  are  very  similar.  It  should  be  noted  that  the 
electrostatic  interaction  energy  decreases  with  increasing  molecule 
size,  especially  in  polar  media  and  when  ionic  species  are  involved. 
Molecule  22  in  Fig.  1  has  a  high  voltage  of  4.52  V  despite  having  an 
IP  value  of  7.33  eV,  and  has  considerably  lower  absolute  solvation 
energy  (-1.28  eV)  than  the  average  AGsol  value  of  -1.85  eV.  The 
low  solvation  energy  is  closely  related  to  the  large  volume  340.9  A3 
for  molecule  22  with  respect  to  an  average  value  of  212.3  A3.  For 
molecule  61  (see  Fig.  2),  however,  the  volume  139.9  A3  is  much 


34 

4.40(7.77) :  4.38 


35 

4.43(7.80) :  4.35 


41 

4.31(7.76):  4.30 


61 

4.23(7.74):  4.16 


Fig.  3.  Eight  molecules  with  IP  >  7.70  eV,  Eox  <  4.50  V.  £0x(calc.):£0x(exp.)  values  reported  along  with  IP  values  in  parentheses. 
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AGS0|  (eV) 


Fig.  4.  Molecular  volume  vs.  AGsoi. 

smaller  than  the  average  value.  Hence,  the  Eox  value  is  only  4.16  V 
despite  a  high  IP  value  of  7.74  eV.  The  molecular  volume  and  AGsoi 
values  are  plotted  in  Fig.  4.  The  correlation  coefficient  is  0.77,  which 
implies  that  the  oxidation  potentials  are  strongly  dependent  on  the 
molecular  size  of  the  organic  molecules  commonly  used  in  LIBs. 

4.  Conclusions 

We  have  calculated  their  IP  and  Eox  values  for  comparison  pur¬ 
poses  with  experimental  values.  The  calculated  potentials  are  in 
close  agreement  with  the  experimental  values,  where  the  RMS  devi¬ 
ation  is  0.08  V  and  the  maximum  deviation  is  0.15  V.  Molecules 
exhibiting  high  Eox  (>4.5  V)  show  one  of  the  following  two  features: 
(1)  IP>7.70eV  or  (2)  IP<7.70eV  with  a  relatively  large  molecule 
size.  Consideration  of  bulk  solvent  effects  is  important  to  describe 
fully  the  experimental  variation  in  Eox‘,  in  particular,  the  electro¬ 
static  attraction  between  solute  and  solvent  is  the  most  crucial 
factor.  Widely  used  PCM  models  reliably  estimate  the  bulk  solvent 
effects  for  the  molecules.  Due  to  its  accuracy  and  reliability,  the 
DFT  calculation  is  recommended  for  virtual  screening  of  the  organic 
molecules  used  as  electrolyte  additives  for  overcharge  protection  of 
LIBs. 


Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  doi:10.1016/j.jpowsour.2008.10.137. 
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